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We show that using the properties of the photon Green's function one can successfuUy describe the 
propagation of arbitrary nonclassical optical radiation through structured materials. In contrast to 
the similar input-output approach, our method is not restricted neither to the spatial homogeneous 
matter nor to the specific direction of light beam incidence or outcoupling. Several quantum states 
of light are addressed in detail, such as Fock, Glauber, and different squeezed light states, and their 
Green functions are given. Furthermore, it is demonstrated how scattering of light at a slab can be 
described with Poynting's energy flow vector, and how experimental setups can be modeled. 

PACS numbers: 42.50.Dv, 42.50.Gy, 42.25.Bs 
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I. INTRODUCTION 

The nonequilibrium photon Green function (GF) 

= (i,(l)ife(2))-(l,(l))(ife(2)). (1) 

is defined in terms of expectation values of correlations 
of the vector potential operator A. Its average A = [ A 
obeys Maxwell's potential equation (wave equation 



□A(r,t) = -^0 [jind(r,i) +jcxt(r,i)] 



(2) 



The photon GF was shown in Ref. [ij to exhibit an exact 
splitting into two components 



Aid - D'^^'P^D-"^ , 



^vac ~ 



ret, — 1 adv, — 1 



'0 



(3a) 
(3b) 
(3c) 



which can clearly be identified as contributions arising 
from electronic generation and recombination processes 
in the medium (P^) in the case of the medium-induced 
contribution -D^^j or, respectively, from the electromag- 
netic state in the free space (vacuum) surrounding the 

medium, represented by the vacuum Green function Dq , 
in the case of the vacuum-induced contribution D^ac- 

Here, elj^* is the transverse dielectric tensor, which de- 
scribes renormalization of a freely evolving vector poten- 
tial 



□Acxt(r,t) = -Aiojcxt(r,t) 



(4) 



due to the presence of the medium to the actual vector 
potential A in linear approximation according to 



ret, — 1 * 



(5) 



If the vacuum GF is given in terms of free plane 

waves, Dva.c can be obtained trivially by replacement of 
the free waves by renormalized waves Aii 



The splitting of [Eq. ([3a|) ] translates directly to the 
energy flux in the system, which is given by the Poynting 
vector, and we can clearly identify the energy flux con- 
tribution caused by the state of light in the free space, 
i.e., its vacuum- induced contribution.^ 

In Secini we will derive and discuss the vacuum Green 
function and the respective vacuum-induced contribution 
to the full Green function for arbitrary states, following 
closely the approach sketched in Ref. [IJ. In particular, 
Fock, Glauber, and different squeezed light states are ad- 
dressed. 

In Sec. IIIH the vacuum- induced energy flow in the slab 
geometry will be developed in more detail to allow even- 
tually for the description of scattering in optical setups 
by Green functions. Modelling of beams splitters and 
photon detection is briefly addressed. 

The advantage of this approach is the exact consider- 
ation of the spatial inhomogeneities inherent to bounded 
media systems. Since the Green functions are globally 
defined, piecewise definition of photon operators is not 
necessary. 



II. VACUUM GREEN FUNCTION FOR 
ARBITRARY QUANTUM STATES 

The vacuum GF consists of the inevitable spon- 
taneous ground state fluctuations of the electromagnetic 
vacuum, which we will denote by -Offsp, and a stimula- 



tion D^g^jjj^ which is given by the preparation of the free 
space to a certain electromagnetic state. 

The complete fluctuations in the free space can be writ- 
ten as 



A<(l:2)-i?0lsp(l-^ 

Note that always Dq^^^^{1,2) 
vacuum spectral function 

Do 



(6) 



^astim(l'2), since the 



^0,sp 



adv 



(7) 



is a unique and universal functionji such that the ^ in- 
dication may be omitted in the following. 
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A. Normal-mode expansion in a mixed state 

We apply in Dq{1, 2) a normal mode expansion for the 
freely evolving vector potential operator— 

Ao(r) = 

Aq V 

I 



where e^q denote the transverse polarization vectors and 
axq (^Aq) are photon annihilation (creation) operators. 
Following the definition ([T|), one obtains the following 
structure: 



Aq A'q' 

+ (C^a.AqA'q' ^^Aq,«(l)^A'q'j (2) + ■ • ■ C.C. . . .)} (9a) 

^o>.,(l> 2) = ^ E E 7=7 {^^n,Vq'Aq FAq.(l)F;,,, ,.(2) + C<,,^,,^ F,^41)F*,^,J2) 

"^^^ Aq A'q' 

+ (a.AqA'q' ^^Aq,,:(l)i^A'q'j (2) + . . . C.C. . . .)} (9b) 



where 



tribution to , 



FAq(l) = exq exp [i(qri - c|q|ti)] 



(10) 



describe classical plane waves with polarization e^q and 
wave vector q. Note that the second term is not the 
complex conjugate of the first to due the index switch. 

The terms cx Cn are normal terms, which depend on 
the difference variables (1 — 2) only and are thus homoge- 
neous. In contrast, the anormalous terms oc Ca depend 
on (1 + 2) and thus pertain to inhomogeneous systems. 
The prefactors C are expectation value differences of cor- 
related and uncorrelated creation and annihilation oper- 
ators in an arbitrary base: 



"-"n,Aq,A'q' 
n, Aq, A'q' 



(^AqftA'. 



("Aq; 



^A'q'. 



{at^ax'q') - (a^J (flA'q') 



Aq/ 



Ca,Aq,A'q' — (flAqOA'q') — (oAq) (oA'q') 



(11a) 
(lib) 
(11c) 



^o>sp,,,(l-2)=i?o<.p,.,(2-l) 

= E^^>..(1)^A..(2), (12) 

Aq 

since, in any (mixed) quantum state described by a sta- 
tistical operator p, 



(flAqaiq,) = Tr{paAqa+q,} 



/ +Tr 



{paA'q'flAq} • (13) 



After separation of Dq^^, the remaining terms in Eq. 
([9|) constitute I?o,stim- They are mutual complex conju- 
gates, and the latter can be written in a more compact 
form by definition of the componentless prefactor 



Application of the commutation relation [a^q, i 



'A'q' J 



'^AqA'q' cxpectation values above (necessary to ob- 

tain normal order in C'^) reveals the ground state con- as 



,AqA'q' = (flA'q'^Aq) " (OAq) (O^-q/ 



(14) 



Do. 



(1,2) 



c 
2iV 



EE 

Aq A'q' 



/qq 



J {C„,AqA'q' FAq,.(l)F;,q, ^.(2) + Ca,AqA'q' ^^Aq.. ( 1 ) ^^A' q' J (2) + . . . C.C. . . .} . (15) 



Using Eq. ([3c]) , the vacuum fluctuations appear renor- cording to^ 
malized due to the presence of a bounded medium ac- 

I?i,(l,2) = i5|(l,2;F^A), (16) 
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i.e., with F replaced by the effective fields 



D. Squeezed light states 



AAq(l) 



.rct,-l 



(l,2)F;,q(2). 



(17) 
(18) 



which describe propagation (i.e., reflection, absorption, 
and transmission) of a classical plane wave in the pres- 
ence of a bounded medium. They are normal mode ex- 
pansions of the vector potential and solutions of Eq. ([5]). 



B. Fock state 

So far we did not make any assumptions on the state 
of the light in the free space. The case of a multi-mode 
number (or Fock) state |{?^}Aq) with the photon popula- 
tion r^;^q = 71^^ ~ fi^^ — 1 was discussed in Ref. |l| already. 
The result from evaluation of the expectation values (fTTj) 
is 



Aq 



+ "!q^A*q,,(l)^Aq,,(2) . (19) 



Uncorrelated expectation values as well as anormalous 
contributions vanish {Dq = Dq). The GF is homoge- 
neous, and its Fourier domain representation is 



TTC 



^ [n%S{u; - cq) + nf^5{u; + cq)] . (20) 



The ground state contribution is contained in the 71^ 
terms and the stimulated part is 



^oW,(l-2)=E^"AqRe [FAq,.(l)i^^*q,,(2)] . (21) 
Aq ^ 



We will now derive and discuss the vacuum GF D< for 



other interesting states. 



C. Coherent state 



For a coherent (Glauber) state with displacement 
P such that^ 



a\f3)^(3\ 
one finds 



|a+ = (/3|/3% ^ <5„,^ , (22) 



Ca = 



(23a) 



A squeezed light state |^, 0) evolves through a finite- 
time interaction, e.g., in a nonlinear crystal, from a given 
state This evolution is described by the unitary op- 
erator S'(^), the squeezing operator: |^, 0) = S'(^) It 
reads in a generalized form for multimode squeezing2i^ 



S{0 = cxp 



EE (^AqA'q'flAqftA'q' 
Aq A'q' 



(24) 



where ^ is a matrix describing the coupling strength be- 
tween different modes. 

The expectation values (jlip can be obtained in a 
straightforward manner using the relations 

S^aXqS = ^ (^/^AqA'q'flA'q' - XqX' q' O-y (25a) 
A'q' 

S^^tci^ = E (/^AqA'q'Ovq. " J^IqA'q'OA'q') , (25b) 
A'q' 

where /i and ly are cosh- and sinh-like power series in 
Obviously, this will lead to an additional mode sum in 
Eq. ([9]), but the general structure remains. 



1. Single-mode squeezing, squeezed vacuum 

The case of single- mode or diagonal squeezing, i.e., 
squeezing of a mode with itself, with ^Aq/2 = 
CAqA'q"^AqA'q'j allows US to quickly obtain the GF of the 
squeezed vacuum. The vacuum GF prefactors for a diag- 
onally squeezed Fock state \{^,n}\ci) are 



Aq'^Aq+I^Aql '^Aq 



"^?~)'5AqA'q' 



CaiO = -AiAqi^Aq(<q + "Aq)^AqA' 



q' 



(26a) 
(26b) 



with the squeezing strength factors fi = cosh |^| and v = 
sinh 1^1 exp(i argf ). The double mode sum in Eq. ([9]) then 
reduces to a single one. We may rewrite 
and find 

C„ =nAq+ |i^AqP + 2|j/AqpnAq (27a) 

Ca = fJ.XqV\q + 2^Aqi^Aq"Aq (27b) 

and decompose the -Do.stim this way into three contribu- 
tions: (i) normal Fock contribution I?^„f [oc n, cf. Eq. 
()19|) ]. (ii) squeezed vacuum contribution [oc alone], 
which remains even for vanishing mode population, 

7?0,sv,.,(l,2) = E 2i^ {kAqP^^Aq,,(l)^A*q,,(2) 



Aq 



-MAq^^Aqi^Aq,^(l)i^AqJ(2) + C.C.} , (28) 
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the squeezed Fock contribution -Do.sf [oc 



and (iii) 

This results was already given in Ref. 1. 



Obviously, 



vanishing of the squeezed contributions -Dq.sv a-nd I?o,sf 
is assured for ^ ^ 0. Note also that these two parts are 
always inhomogeneous, since they contain a Ca term. 



2. Squeezed vacuum in a mixed state 

Since any mixed state may be expanded into Fock 
states, 



Tr{p$} = ^(n,|$|n,) 



(29) 



it is easy to see that the squeezed vacuum GF, Eq. 
will appear in any single-mode squeezed (mixed) state. 
Evaluating this Fock state expansion for general multi- 
mode squeezing reveals a squeezed vacuum GF as well, 
with the same structure (but different definition of fijv). 
In conclusion, any squeezed light GF exhibits a squeezed 
vacuum contribution of the form (|28p . 

Note again that Z?o,sv is defined here such that it is 
free from contributions from the ground state, Do.sp- 



3. Single-mode squeezed coherent light 

For a diagonally squeezed coherent state 
correlated and uncorrelated expectation values largely 
compensate in ZJo.stim, and only the squeezed vacuum 
remains. 



C'a, AqA'q' — ^MAqi^Aq<5AqA'q' 



4- Two-mode squeezing 



(30a) 
(30b) 



In a widely used setup for the generation of squeezed or 
entangled light, the optical parametric oscillator (OPO), 
photons are emitted from a nonlinear crystal in pairs with 
correlated moments qo ± q around the reference momen- 
tum qo. It may be regarded as an outer parameter of 
the systern and q as a variable. Usually, the squeezing 
strength ^ is assumed to be symmetrical around qg. It 
vanishes for q = and has a finite bandwidth 
for q — > oo). 

The squeezing operator can be obtained from Eq. ((24|) 
by diagonalization:— 



Sii) = exp 



(CAq^A,qo+q'JA,qo-q 

Aq 



-^Aqatqo+q«A,qo-q) (^1) 



For this case, it is convenient to express the coefficients 
fi, v hy functions 



MA(|qo - q|) = MAq , I'Adqo - q|) = i>Aq, 

such that the relations (1251) read now 



aXqS — Z-iAqflAq - t^AqaA,2qo 



^ "-Xq^ — /^AqflAq ^ i^AqaA,2qo-q ■ 



(32) 

(33a) 
(33b) 



In a Fock state, the expectation values are [cf. Eq. ((26)) ] 



C|(0 = (/^Aq^Aq + I^AqP?lJ,2qo-q)W'q' 
C'aii) = -AAqJ>Aq(nA,2qo-q + ^iA,2qo -q ) W'q' 

or, in the compact notation for -Do.stim alone: 

C"(|) = (It^AqP + ?^Aq + It'AqP^Aq 

+ pAqP"-A,2qo-q)'5AqA'q' 
Ca(0 = ("MAql^Aq ^ 2/iAq?'Aq?^A,2qo -q) SxqX'q' 



(34a) 
(34b) 



(35a) 
(35b) 



Again, the GF decomposes into three components. Com- 
paring them to the diagonally squeezed Fock state, one 
finds that, except for /i ^ /i, only the Z^o.sf component 
differs. 



E. Propagation of squeezed light 

The propagated squeezed light correlations (and their 
propagated components) may be obtained as before ac- 
cording to Eq. (fT6| . Their complex prefactors C„ , Ca and 
inhomogeneities introduced by squeezing remain. 



III. ENERGY FLOW WITH NONCLASSICAL 
LIGHT 

We will now consider energy flow and light scattering 
with the help of Poynting's theorem;^ dU/dt -\- divS = 
— jE, following closely the lines of Ref. [H but extending 
and refining the approach presented there. 

Up to now, the results presented in this work are valid 
for arbitrarily inhomogeneous and instationary systems. 
For an easier discussion and better comparison with for- 
mer work, we will now regard a system in slab geome- 
try, i.e., TE-polarized light propagating along the x axis 
through a linear medium slab of thickness L which is in- 
finitely extended in the y-z plane, isotropic and steadily 
excited. Incoming light may be instationary; its polar- 
ization is chosen along the z axis. 

The vector potential now takes the form 

ATE,q(r,t) = e^exp [iq||r|| - icqt] Aq{x) , (36) 

where qy and are the in-plane components of q and r, 
and q± and x the corresponding x direction components. 
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In contrast to, e.g., Ref. |l|, we will need to properly The splitting of the GFs translates directly to the in- 

consider the coherent contribution ((A) ^ 0) to the en- coherent component of the energy flow. For the spon- 

ergy flow. In slab geometry, energy flow is only possible taneous contribution to S^, i.e., due to D^,.^ [Eq. (fn])]. 

in the x direction, and the corresponding Poynting vector obtains^ 
component reads after insertion of the photon GFs 

5,(1) = 5,,ich(l) + ^.,coh(l) (37a) 



2-»l 



5.,co.(l)--^|^^(A.(l))(A.(2) 



(37b) 
(37c) 



2^1 



The contribution from Do.stim in its general form [Eq. p5|) ]. 

he? Fq ( d r 

5'x,stim(l) = — lJ-^\\—i^^ ( C„ylq(a;)— A*,(a;)exp i(q|| - qj, )r|| - ic(g - q')* 



q.q' ' 



~CQAq(a:)— Aq.(a:)exp i(q|| +q||)r|| -\c{q^cl)t 



(39) 



exhibits stationary terms only for q' = q. In the case of squeezed vacuum follows from C„ = |i^AqP<5AqA'q' [Eq. (|28l) ] 

5x,sv(l) = 5Z ( l^ql'l"i^q(^)^^q(2^) +Aiq|i'q|Ini |Aq(x)— Aq(a:) cxp [2i(q||r|| -c(7t) + iargO] I j • (40) 

Its homogeneous terms are stationary, the inhomogeneous terms oszillate with double frequency. 

A similar analysis shows that Sx,coh is formally identical to S'^.stim, so that both can be treated the same way in 
the following. 



'S'2;,coh — '5'a;_stim[C'n ^ (o-q) (tlq') : Ca ^ (o-q) (^.q')] 



(41) 



A. Scattering of incident light 

The theory so far can represent any field in the free space, since the vacuum Green function is a sum over all 
wave- vectors q, and an arbitrary state |0) is allowed. We may now address a specific situation by assuming a specific 
structure of the vector potential. 

The following vector potential represents light incident onto both the left and the right side of the slab {symmetric 
setup, with complex amplitudes ql and gj^), which is reflected and transmitted to the respective other side: 



Aq(x) 



gj^e-'i^^ + [gnr^ + gLt^)e<^^^ x > +L/2 



(42) 



where rq, tq are complex reflection and transmission coefficients, respectively. The vector potential Aq for coordinates 
inside the slab is not used in the following and does not have to be known. This trivially covers the case of multiple 
incident light beams of wave vector q, which superimpose to give a new complex amplitude gL/R- One would set 
gu = to obtain the rcficction and transmission of light incident from the left alone {asymmetrie setup). 
The energy flow density through the medium surface (area F) is 



AS* 



div S dV 



1 

F 



Sdf^S4L/2)-Sxi-L/2). 



(43) 



Here, = Sx{x> ^/2) describes the energy flux emitted from the right surface into the free space in parallel to the 
x axis, while = Sx{x < — V^) is the left surface counterpart emitted into the negative x direction. 
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For the calculation of S^/^ according to Eq. p9|) . the following expressions have to be evaluated: 



= Tig± (|.9L/flP - \9LlRr^f - Igfl/iiqP + 2ilm {I^l/aPK + iq) exp2ia;}) (44) 



X> + -f-/2 



The Im term in the former equation does not contribute to S^l^ if ImC„ = 0, as, e.g., in (squeezed) Fock states. 
Let us first analyze the spontaneous contribution (where always \gh\ = |.9_r| = !)• We obtain 

q 

The balance of these energy flows is, as it should, cancelling out in the symmetric setup {AS — 0). To analyze the 
scattering, we set gn = and readily reproduce the result of Ref. Ground state fluctuations are transmitted and 
reflected just as classical light, with (l — |rp — = a, where a is the classical absorptivitypii^ with amplitude 1/2. 

The -Do.stim contribution to S*^/^ can be obtained by straightforward calculation. We concentrate on cases with 
Cn.a oc ^AqA'q' and Im C„ = for simplicity. With the help of the asymmetric setup scheme, we find that the energy 
fiow of the incident light in general reads 

SstL = — E'?-Ll5iP (C» - ICal cos[2q±x-2cqt + ipq]) (p^ = 2arggi + argC^ , (47) 

q 

while the energy flow after transmission to the right is (for reflected light accordingly) 

^stL = — E^^l^^l^l^il^ ^ |Ca|cos[2(jf_La; - 2cqt + Lpq\) = 2 arg + arg Ca + 2argiq. (48) 

q 

i.e., the entire energy flow is damped by |iqp, and a phase delay of 2 arg tg in the instationary part is accumulated. 



B. Beam splitters and photon detectors 

A beam splitter, e.g., could be described following this 
scheme as a slab onto which light is incident from both 
the left and the right under an angle of 45 degrees. Then, 
the linearly superimposed transmitted and reflected en- 
ergy flows on each side are the output channels. Any 
absorbing or (spatially) dispersive behavior of the beam 
splitter is fully accounted for by the present approach. 

A photodetector measuring an electromagnetic mode 
converts photons in photoelectrons, hence giving rise to 
an electric current, called photocurrent i. Since the 
Poynting vector gives the photon energy per second and 
square meter, it is natural to assume that the mean 
value of this photocurrent recorded during the (small) 
time interval t,t + At is proportional to the normally 
ordered energy flow (: S :) of the incident photons, i.e., 
{i{t,At)) — rjAt {: S{t) :), where 77 is the quantum effi- 
ciency factor of the detector. Normal ordering eliminates 
the ground state fiuctuation contribution [cf. Eq. (fTS]) ]. 
which is usually not measured by detectors. It is equiva- 
lent to a restriction to the stimulated and coherent con- 
tributions in our approach. Hence, 

{■.S:)= 5,tim + 5coh , (49) 



is the quantity that determines the photodetection mea- 
surement outcome. 



IV. CONCLUSION 

In this article, we show in detail how the Green func- 
tion approach to the scattering of squeezed light pre- 
sented in Ref. [l] can be generalized to light in arbitrary 
quantum states. 

Green functions for several quantum states in the free 
space are derived. They can generally be separated into 
a ground state fluctuation and a stimulated contribution. 
In particular, squeezed light is addressed and the cases 
of general multimode squeezing, single-mode squeezing, 
and two-mode squeezing are discussed. All squeezed light 
states are shown to result in the same squeezed vacuum 
contribution as a part of the stimulated contribution to 
the vacuum GF. 

Vacuum GFs for arbitrary quantum states are con- 
structed after normal-mode expansion of the vector po- 
tential. They are represented by normal and anor- 
malous (inhomogeneous) terms with appropriate prefac- 
tors Cn , Ca given by photon operator expectation values 
in the considered quantum state [Eq. ©]. 
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The normal modes propagate as classical waves in any 
(bounded) media system. Thus, these states are scat- 
tered no other than classical light. This conclusion is 
possible thanks to the exact splitting of the Green func- 
tion into medium- and vacuum-induced contributions^i 

The same property allows to consider independently 
the electromagnetic energy flow contribution caused by 
light originating from the free space. Consequently, we 
develop the Poynting vector for stimulated light in its 
general normal-mode representation [Eqs. ([57|) - (|¥l]l ]. and 
then discuss the description of light propagation and scat- 
tering through a slab on this basis. The advantage of 
this approach is the exact consideration of spatial inho- 
mogeneities in the system and its validity for oblique in- 
cidence. Other than in the input-output formalism, there 
is no need for a spatial decomposition of the photon op- 
erators into incoming and outcoming, since the GFs are 
global functions. 

After propagation through a medium, the energy flow 



is damped by |ip and a phase delay of 2argi is accu- 
mulated [Eq. (pS)) . for reflected hght accordingly]. This 
simple result is proven here for arbitrary light states and 
arbitrarily absorbing or dispersive media. 

Finally, it is shown how to describe some elements in 
an experimental setup with the present approach, for 
which the beam splitter serves as an example, and how 
the Poynting vector as a natural quantity for the energy 
flow relates to the outcome of measurements with pho- 
todetectors. 
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